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A  B  S  TRACT 


The  transient  and  steady-state  behaviour  of  the  M/M/r 
queueing  process  with  bulk  arrivals  is  analysed.  The  transient 
behaviour  is  treated  in  terms  of  Laplace  transforms,  and  steady- 
state  behaviour  -  in  terms  of  generating  functions  and  probabili¬ 
ties.  The  influence  of  the  bulk-size  variance  on  the  expected 
queue  size  is  discussed  at  some  length. 


MULTIPLE-SERVERS  QUEUE  WITH  BULK  ARRIVALS 


D.2.  Mittwoch  and  B.  Avi-Itzhak 


INTRODUCTION 


In  many  queueing  situations,  customers  arrive  in  bulks  and  are  served 
individually.  Such  is  the  case,  for  example,  in  an  airport,  where  passengers 
arrive  as  a  single  group  in  a  plane  but  are  served  individually  at  the  passport- 
control  and  customes  counters.  The  present  paper  deals  with  a  system  where  bulks 
with  randomly  distributed  size  arrive,  in  a  stationary  Poisson  stream,  at  a  single 
queue  attended  by  r  servers,  the  queue  discipline  being  first-come-first-served. 
Members  of  a  newly-arrived  bulk  join  the  end  of  the  queue  in  a  random  order. 
Service  times  are  assumed  to  be  mutually  independent,  with  an  identical  negative 
exponential  distribution.  The  transient  and  steady-state  behaviour  of  the  system 
is  analysed,  and  some  closed-form  results  concerning  queue  sizes  are  derived. 

A  similar  bulk-arrivals  model,  in  which  service  is  provided  in  batches, 
was  derived  by  Loris-Tegham  [1]  who  formulated  the  Kolmogorov  equations  and 
obtained  solutions  for  special  cases.  The  case  r  -  1,  a  single-server  system, 
was  studied  extensively  by  Gaver  [2]  using  the  imbedded  Markov  chain  method. 

Here  the  solution  is  approached  through  Kolmogorov's  forward  differential- 
difference  equations,  using  generating  functions  and  their  Laplace  transforms 
In  the  steady-state  situation,  particular  attention  is  given  to  the  effect  of  the 
mean  and  distribution  of  the  bulk  size  on  the  average  queue  size. 


2 


MATHEMATICAL  MODEL 


a.  Transient  behaviour 

The  bulk  arrival  rate  is  denoted  by  A  and  bulk  size  is  assumed  to  be 
arbitrarily  distributed  as  a  random  variable  N,  possessing  a  probability  mass 
function  fN(n),  n«l,2,...  .  Tne  service  rate  for  each  of  the  r  servers  is 
denoted  by  p. 


Let  Pji^  denot«  the  probability  of  there  being  i  customers  in  the 
system  at  time  t,  given  that  at  t«0  their  number  equals  j.  For  simplicity, 
we  omit  the  subscript  j  and  use  P^t)  instead  of  Pji(t).  Kolmogorov's 
forward  differential -difference  equations  take  the  following  form: 


3t  pi(t)  *  Xm5o  fN(i-ffl)Pm(t)  *  Min{i*l,r)uPi+1(t) 
-  (A+min{i,r}p)PiCt)  ,  i»l,2,,..  , 

and 

at  po(t)  -  -  xp0(t)  4  »piW  • 

Note  that  f^(n)  ■  0  for  n<l. 

We  now  define  two  generating  functions: 


U) 

(2) 


G(z,t)  «  i|0Z1Pi(t)  , 

■  nfozVn>  • 

The  Laplace  transforms  of  G(z,t)  and 
G*(z.s)  -  |e"StdG(z,t)  , 


1*1  f  1  <• 

1*1  f  1  • 

P^(t)  are  defined  as 
Re(s)>0  j, 


Pi(s)  *  Ie’StdPi^5  > 


(3) 

(4) 

(5) 


Re(s)>0  . 


(6) 
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The  Laplace  transform  G*(z,s)  is  obtained  from  Eqs.  (1)  and  (2)  by  taking 
first  the  generating  functions  and  then  referring  to  relations  (S)  and  (6): 

z^+d-Du^r-iJ^PMs) 

G  ^z,s)  *  sz+ta(l-GN(z))  +  (z-l)ru  ’ '  ^ 

where  j  is  the  number  in  the  system  at  time  t«0. 

In  Eq.  (7)  the  Laplace  transforms  P*(s),  i«0, 1  ...,r-l  are  unknown. 

The  following  method  is  suggested  for  determining  these  functions:  define  a 
generating  function 

Q(z.t)  -  ^£p*iPi(t)  ,  |z|  <  1  .  (8) 

From  Eqs.  (1)  and  (2),  we  obtain 


§£  Qd.t)  -  w(l-z)  fj  Q(z,t)  - 


AijlzlJofN(1’B)Pm(t)  *  xQ(*»t)  +  (r-l^zr‘2Pr.j(t) 


(9) 


This  nonhomogeneous  partial  linear  differential  equation  can  be  solved  for 


in  terms  of 

p^t),  1*0,1, 

•  •  • »r-2. 

The  solution  is  obtained  in  the  form 

P(u(z,t,Q), 

v(z,t.Q))  -  0, 

where 

and 

u(z,t,Q)  - 

C1  * 

Cj  *  constant  , 

satisfy  the 

v(z,t,Q)  ■ 

equations 

C2  ' 

Cj  ■  constant  , 

dt 

T  * 

dz 

-  *«»»*> _ 

1  W(Z‘1)  ^QC*.t:)-xJ|Jzi^|AfNCi-in)Pi|i(t)-(r-l)uzr"2Pr  jCt) 


Q(z,t) 


(10) 


The  functional  form  cf  ff  is  obtainable  from  the  initial  condition  of  the  process 
(For  details  see  Sneddv ..  [3] ;  the  same  method  is  used  by  Saaty  [4]  in  solving  the 
M/M/r  process). 


After  some  rather  lengthy  manipulations,  the  desired  solution  is  obtained 
in  the  following  form: 
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where 


Q(*,t) 


’  Ifl  l&  {V(tHV( l-.).-,,(t-u))1£I|(i^)P.(u)du 

♦  Q(*,0)*(r-l)/teX(t*u)(l-(l-i)e‘uCt*u:))r'2P  ,(u)du 
o  r- 1 


(ii) 


'0 

Q(i,0)  - 

,  zJ 


if 


j  >r-2 

otherwise  . 


Taking  the  Laplace  transfora  of  Eq.  (11),  we  have  (note  that  the  right- 
hand  side  of  Eq.  (11)  is  given  in  foras  of  convolutions): 


fu(i-»)6(*-l)i"' 


rSV  p*fsl  -  rt2  iI1  i  NU'aJ  [wHI 
i-0  i  i-1  a*0  w«0  s+(i-w)p+X 


p;(s) 


r-2  (r;2)(*-J)r‘2*W 
*  QU.OJ^r-l)^  -4-;(r-z-w)M;\  p;.iW 


(12) 


Equation  (12)  is  an  identity  of  two  polynomials  of  degree  r-2  in  z.  The 
coefficients  of  equal  powers  of  z  on  both  sides  aust  be  identical,  hence,  (r-1) 
linear  equations  are  obtainable,  relating  P«(s),  i-0,1, . . . ,r-l;  these  are 
conveniently  derived  by  repeated  differentiation  with  respect  to  z  it  s  -  0, 
and  yield 


n, ,n-i. 


"!<■’  ■  sr.  ili  X  &-•>'  1  >:<•> 


r-i-2. 


i-0,1 . r-2. 


(13) 


where 


T .[1  ^  i • J 

[  0  otherwise. 

Note  that  in  Eq.  (13),  and  throughout  this  paper,  •  sua  in  which  the  upper  limit 
is  less  than  the  lower  Unit  is  defined  as  equal  to  zero. 

We  now  have  r-1  linear  equations  and  r  unknowns,  namely  PJ (s) , 
i-0,1, . . . »r-l.  The  missing  relation  is  obtained  froa  Eq.  (7).  G*(z,s)  is 
analytic  with  respect  to  z  in  |z|<  1.  The  denominator  of  Eq.  (7)  has  exactly 
one  zero  in  |z|<  1,  which  must  also  be  a  zero  of  the  numerator  yielding  the  miss¬ 
ing  linear  equation  in  P|(s),  i-0,1,. ...  r-1.  To  prove  that  our  dsnoeiinator  has 


s 


one  zero  as 

above,  let 

g(z) 

■  •  AzGn(z) 

-  tv 

and 

f(0 

■  (A  ♦  ru  ♦ 

s)z 

For  |z |  ■  1 

we  have 

1*00  1 

>  lg(*)l  • 

It  follows  then  from  Rouche's  Theorem  that  f(z)  ♦  g(z)  (which  is  our  denomina- 
tor)  and  f(z)  have  the  saae  number  of  zeros  in  z<  1,  and  since  f(z)  has 
exactly  one  zero  in  |z|<  1  (at  z«0) ,  the  denominator  of  (7)  has  exactly  one 
zero  in  |z|<  1.  On  |z|  ■  1  the  denominator  cannot  equal  zero,  since 

|f(z)  *  g(z) |  >  |f(z)  |  -  |g(z)  |  >  0  , 
and  the  proof  is  complete. 

The  single  zero  in  question  is  real  for  s  ■  Re(s) .  If  we  denote  it  by  zq 
the  additional  linear  equation  takes  the  form 

<vi}  ♦ iiol - ° 


Clearly,  only  numerical  solutions  are  possible  for  the  general  case,  since  zq 
is  unobtainable  in  closed  form.  For  bulk  sizes  not  exceeding  3,  zq  can  be  deter 
mined  in  closed  form  as  a  solution  of  a  polynomial  whose  degree  is  the  maximum 
bulk  size  plus  one. 


l.ie  form  of  P|(s),  i>r,  is  readily  obtained  from  Eq.  (7) 


pi<*>  ■  [4H£ttt2 


i-j 

otherwise, 


where  n 


if 


(15) 


A(z,w)  -  ^  (sz4Az(l-GN(z))  +  (z-l)ru)'1  , 
9zW  " 

and  A(z,w)  ■  0  if  w<0. 
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b.  Steady-State  behaviour 


The  steady-state  generating  function  for  the  number  in  the  system,  denoted 
by  G(z),  is  obtained  as  follows: 

G(z)  -  |fi0  sC*(*.s)  ■  nTa-^-n.a-Orp  •  0*) 

where  P.  ■  lia  P, (t) . 

The  values  of  P^,  i«0,l, . . . ,r-l,  are  obtainable  as  the  solution  of  the 

following  set  of  r  linear  equations: 

■  g  v*~>v  a  a  $<-»'  — 

(rt‘i 


<r-1)»Pr-l<r;!)rti;2(-l)'  -rht  .  l-O.l.! . r-2  ,  (17) 


and 


"Ufrr-DPi  “  r“-lEW- 


(18) 


The  (r-1)  equations  given  by  (17)  were  obtained  as  the  liait  of  equations  (13), 
and  the  r-th  equation,  (18),  is  obtained  by  setting  z«l  in  equation  (16).  It 
can  be  shown  that  the  steady-state  condition  is 


p  .  iM<i  . 

ru 

For  example,  for  r  ■  3,  the  solution  of  (17)  and  (18)  is: 
P  ■  - 


(19) 


3  -  rp 

ir~r t 


3  *  2  5  ♦  fl  <5  ♦  1  -  fN(l» 


1 


x.  a  x 


(201 


P1  ■  5  po  •  *nd  p2  ■  2  ((1  +  p  5  - ;  • 

The  expected  maiber  of  customers  in  the  systea,  denoted  by  L,  is  readily 
obtained  from  Eq.  (16)  as 

2 

1  '  S  ((| $j-  ♦  ♦  2i(l-p)) (r-i)Pi  (21) 
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«nd  the  variance  of  the  number  in  the  system,  denoted  by  V,  is  obtained  as: 


..  .  nl  *  111  ■  D  •  3i(i-i)ti-»)  » 


.’(f®  •  »2 


3r(l  -  p) 


2r(l  -  p) 


T  } 


(r-i)Pi  ♦  L(l-L) 


(22) 


In  the  special  case  where  N  is  a  constant,  say  N  -  C,  Pi#  i»0,l,...,r-l, 
is  obtainable  in  closed  fora  as  follows: 


and 


- 1  r(1  ~  V  . 

r  *  nSl  (-nT  .Bo  ♦  ■» 


o  isl 

n  .So 


(-  ♦ 
p 


i-1,2 . r-1 


(23) 


Substitution  in  (21)  and  (22)  yields  closed-fora  expressions  for  L  and  V. 

Figure  1  gives  the  value  of  L/r  as  function  of  the  expected  bulk  site, 
E(N),  for  different  values  of  r  with  p  kept  constant  at  0.9.  The  solid  lines 
represent  the  case  of  constant  bulk  size  and  the  dashed  lines  that  of  geometrical¬ 
ly-distributed  bulk  size.  The  figure  shows  that  L  increases  approximately 
linearly  with  E(N) .  Also,  for  the  geometrical  distribution  L  is  greater  com¬ 
pared  to  the  constant  bulk  size  case. 

Intuitively,  one  would  expect  L  to  be  smallest  when  the  bulk  size  is 
constant,  thav.  is  f^(E(N))  ■  1.  This  conjecture  is  also  supported  by  Figure  1, 
where  L  for  the  constant  bulk  size  was  found  to  be  smaller  than  for  its  geometric¬ 
ally  distributed  counterpart  with  identical  Bean.  Here,  only  one  sufficient  con¬ 
dition  will  be  proved,  a  more  general  proof  being  too  complicated. 


Theorem:  Let  A  and  B  be  two  queueing  systems  of  the  type  described  in  this 
paper,  identical  in  all  respects  except  for  the  bulk-size  distribution.  If  the 
bulk  sizes  in  A  and  B  possess  probability  mass  functions  *WA(-)  «nd  (•) 

respectively,  such  that  ”(N^)  -  E(Ng)  -  C  >  r-1,  and  fN  (i)  -  0  for  1-0,1 . r-2, 

then  La>Li,  if  «nd  only  if  V(NA)>V(NB). 
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Proof.  Let  PJA)  and  P^B),  i-0,1,2,...  denote  the  steady-state  probabilities 
of  there  being  i  customers  in  system  A  and  B,  respectively.  From  Eqs.  (17)  and 
(18)  it  follows  that: 

(A)  *  »  i«0,l,...,r-l  .  (24) 

Substitution  of  Eq.  (24)  in  (21)  yields: 

V(N  )  -  V(N.) 

LB  '  LA  '  “JTcTTt^r-  Jo  Cr-DPJA)  .  (25) 

From  Eq.  (25)  it  is  seen  that  (Lfi-  LA)  is  an  increasing  linear  function  of 

(V(Nb)  -  V(Na)).  For  the  special  case  where  the  bulk  sire  in  system  A  is  constant, 
we  have: 


V(Nb) 

lb  -  la  *  gen:p)!t 


$  tr-OPjA)  > 


(26) 


EM 

Figure  /.  Expected  number  of  customers  in  the  system 
per  server  as  a  function  of  the  expected 
Au/A  stye ,  (?*o.9). 
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